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Aims. Fundamental properties of the extinction curve, like the slope in the rest-frame ultraviolet (TJV) and the presence/absence of 
a broad absorption excess centred at 2175 A (the UV bump), are investigated for a sample of 108 massive, star-forming galaxies at 
1 < z < 2.5, selected from the FDF Spectroscopic Survey, the K20 survey, and the GDDS. 

Methods. These characteristics are constrained from a parametric description of the UV spectral energy distribution (SED) of a galaxy, 
as enforced by combined stellar population and radiative transfer models for different geometries, dust/stars configurations and dust 
properties. 

Results. In at least one third of the sample, there is a robust evidence for extinction curves with at least a moderate UV bump. The 
presence of the carriers of the UV bump is more evident in galaxies with UV SEDs suffering from heavy reddening. We interpret 
these results as follows. The sample objects possess different mixtures of dust grains and molecules producing extinction curves in 
between the average ones of the Small and Large Magellanic Cloud, where the UV bump is absent or modest, respectively. Most 
of the dust embeds the UV-emitting stellar populations or is distributed out of the galaxy mid-plane. Alternatively, even dust with 
a pronounced UV bump, as for the average Milky-Way extinction curve, can be present and distributed in the galaxy mid-plane. In 
this case, variations of the continuum scattering albedo with wavelength or an age-dependent extinction are not sufficient to explain 
the previous trend with reddening. Hence, additional extraplanar dust has to be invoked. The data suggest that the carriers of the UV 
bump are associated with intermediate-age stellar populations, while they survive in the harshest UV-radiation fields owing to dust 
self-shielding. 

Conclusions. The existence of different extinction curves implies that different patterns of evolution and reprocessing of dust exist at 
high redshift. Ignoring this may produce a non-negligible uncertainty on the star-formation rate estimated from the rest-frame UV. 
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1. Introduction 

Extinction occurs whenever electromagnetic radiation propa- 
gates through a medium containing a mixture of dust grains and 
molecules. The dependence of absorption and scattering by dust 
from the wavelength of a photon is described by the extinction 
curve, which depends on the physical and chemical properties 
of dust grains and molecules (Whittet 2003). On a galaxy scale, 
the propagation of electromagnetic radiation through a dusty 
medium is described by the attenuation function. This is the 
combination of the extinction curve with the geometry of the 
system, in which a substantial fraction of the scattered light is 
returned to the observer's line of sight. 

The direct determination of the extinction curve from obser- 
vations is limited to our own Galaxy, the Magellanic Clouds, 
and M 31 (e.g., Fitzpatrick 2004, Clayton 2004). The most strik- 
ing differences among the typical extinction curves of the Milky 
Way (MW) and the Large and Small Magellanic Clouds (LMC 
and SMC, respectively) occur in the UV domain, where small 
dust grains and molecules absorb and scatter with the highest 
efficiency. From the MW to the SMC, the presence of a broad 
absorption excess centred at 2175 A (the UV bump, see Witt & 
Lillie |1973] and references therein) vanishes (almost) completely, 



Send offprint requests to: S. Noll 



while the slope of the extinction curve in the far-UV becomes 
steeper 

Multiple carriers (i.e., organic carbon and amorphous sili- 
cates) were suggested to explain the enigmatic, invariant cen- 
tral wavelength and variable bandwidth of the 2175 A feature 
(Bradley et al. 2005 and references therein). The ingredients of 
interstellar dust that originate the UV bump and the far-UV rise 
are in broadly similar proportions in the MW and LMC, but dust 
seems to be more radically different in the SMC (Whittet l2003l l. 
Differences of the same relevance exist within our own Galaxy 
for different sight lines and towards different environments, from 
dense molecular clouds to the diffuse interstellar medium (ISM; 
see Fitzpatrick 2004). They must reflect the extreme sensitivity 
of the small-sized dust components to the local chemical enrich- 
ment and energy budget (in terms of radiation field and shocks), 
as well as their selective removal from the size distribution ow- 
ing to a number of physical processes (see Whittet 120031 ). 

Nearby starburst galaxies have an extinction curve that 
seems to lack a 2175 A bump, as the SMC curve, and to have 
a steep far-UV rise, intermediate between the MW and SMC 
curves (Gordon et al. I1997I I. This extinction curve of starbursts 
was constrained from the modelling of the observed UV/optical 
broad-band colours. An alternative interpretation of the empir- 
ical attenuation law in starbursts (Calzetti et al. |1994| l invokes 
a turbulent interstellar medium (Fischera et al. 120031 ). However, 
it adopts a foreground screen for the distribution of MW dust 
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that is reminiscent of the dust/stars configuration of the Gordon 
et al. (1 19971 1 model that reproduces the so-called "Calzetti law". 
This configuration is supported by models of large-scale galactic 
winds driven by momentum deposition (Murray et al. |2005l l. 

SMC-like dust seems to be present in Lyman-break galaxies 
at z ~ 3 as well (Vijh et al. 2003 ). This result is complemented 
by the finding of Noll & Pierini (2005. INP()5] l. based on spec- 
troscopy instead of broad-band photometry. For a sample of 34 
massive, UV-luminous galaxies at 2 < z < 2.5, INP05] find that 
the majority of the objects with strongly reddened, rest-frame 
UV spectra have an extinction curve similar to the LMC curve. 
Conversely, the objects with the least reddened, rest-frame UV 
spectra seem to host SMC-like dust. This result strengthens the 
evidence for a difference in the properties of the dusty ISM at in- 
termediate/high redshifts arising from previous, sometimes con- 
tradictory, results in the literature (e.g., Malhotra 1997; Pitman 
et al.l2050l Maiolino et al. 2001 , Vernet et al 2001 , Hopkins et 
al.'2004'i Savaglio & Fall |2004i Wang et al. , 20041 Wild & Hewett 
[2005; Yorket al. .2006) . 



Here we extend the study of |NP05l to a heterogeneous sample 
of massive, UV-luminous galaxies at 1 < z < 2.5. This sample 
is three times larger than the INP05l one and probes the new red- 
shift range 1 < z < 2. Furthermore, the available ground-based 
spectroscopy allows us to date the oldest stellar populations at 
least in objects at 1 < z < 1.5. We can also characterise the mor- 
phology of a galaxy across the whole redshift range in a reliable 
way from existing imaging with the ACS camera on board the 
Hubble Space Telescope (HST). 

Throughout this paper Hq = 70kms ' Mpc"', Oa = 0.7, 
and Q.M = 0.3 are adopted. Photometry is given in the Vega 
magnitude system. 



2. The spectroscopic sample 

Inferring the presence of the broad dust absorption feature cen- 
tred at 2175 A either directly (from individual and stacked spec- 
tra) or indirectly (with the method established in INP05I see 
Sect. 13.11 1 from galaxy spectra taken at optical telescopes lim- 
its the redshift range to 1 < z < 2.5. The need for rest-frame 
UV spectra with a good continuum definition limits the sample 
to bright (and massive), star-forming galaxies. Finally, the need 
for a representative sample of such galaxies within this redshift 
range implies selection from different spectroscopic surveys. 

For the total sample of 108 massive, star-forming galaxies 
at 1 < z < 2.5 selected as described hereafter. Figs. [T] and 
|2] show the distribution of the R magnitude (observed frame) 
as a function of redshift and the distribution in the R - Ks-K^ 
colour-magnitude plane (observed frame), respectively. It is ev- 
ident that our sample contains mostly objects with 23 < /? < 24 
and 2 < R - Ks < 5 whatever the redshift, R - Kf^ turning bluer 
towards fainter magnitudes. In general, the constraint on the 
S/N of the optical spectra implies that galaxies at 1.5 < z < 2.5 
are more luminous in the rest-frame UV than those at 1 < z < 1 .5 
(cf. Fig. [8]l, on average. In terms of star-formation activity and 
stellar mass, the total sample is less heterogeneous than what 
may be expected from the different selection criteria of the three 
surveys of origin (cf. Sect. |4.4t . Table[T]lists basic properties of 
the sample galaxies. Hereafter we illustrate the different selec- 
tion criteria and properties of the three subsamples. 
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Fig. 1. R magnitude versus redshift for the total sample of 
108 actively star-forming galaxies at 1 < z < 2.5 selected 
from the FDF Spectroscopic Survey (lozenges and crosses), the 
K20 Survey (circles and plus signs), and the GDDS (asterisks). 
Crosses, plus signs, and asterisks mark galaxies without a deter- 
mination of the UV continuum slope. In order to illustrate se- 
lection biases, the curves show lines of constant luminosity for 
typical 1 < z < 1 .5 FDF (solid line) and K20 (dashed line) spec- 
tra. Both curves intersect atR- 23.5 and z - 1.5. 
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Fig. 2. R - Ks versus magnitude for the total sample under 
investigation. Symbols are the same as in Fig.[T] 



2. 1 . The FDF subsample 

In addition to the 34 galaxies at 2 < z < 2.5 with R < 24.7 inves- 
tigated bv lNP05l we select 32 objects at 1 < z < 2 with R < 24 
from the same /-limited FORS Deep Field (FDF) Spectroscopic 
Survey (Noll et al.'2004). Only galaxies with a dominant young 
stellar population (types III to V, see Noll et al. 2004) are consid- 
ered. Basic aspects of the data reduction are illustrated in NoU et 
al. (E004i) and lNPOSl 

Ground-based optical/near-IR photometry in nine filters is 
available from Heidt et al. ([20031) and Gabasch et al. (2004^, re- 
spectively. The FDF was also imaged in the broad-band F814W 
filter with the HST ACS camera. A lOcr limit of 26 mag was 
reached with four WFC pointings of 40 min-exposure each. 
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The data reduction was performed with the standard CALAC^ 
pipehne, and the combined final mosaic was produced with the 
MultiDrizzle package (Mutchler et al. ,2002j . 

2.2. The K20 subsample 

We select 34 star-forming galaxies at 1 < z < 2.3 from the spec- 
troscopic catalogue of the K20 Survey in the Chandra Deep Field 
South j;CDFS) and a field aroun d the quasar 0055-2659 (Cimatti 
et al. 1X002; MignoH et al. |2005] l. Besides having K, < 20, all the 
selected objects but two have R < 24. This subsample comprises 
only K20 objects with a secure redshift and a suitable spec- 
tral coverage across 2175 A (basically from 1900 to 2450 A). 
Galaxies with spectral energy distributions (SEDs) dominated by 
old stellar population (class 1 and 1.5, see Mignoli et al. |2005l l 
are rejected. 

In addition to the basic reduction described in Mignoli et al. 
( 1200 5'), the K20 spectra have been corrected for slit losses using 
the R magnitude of the individual galaxies. Furthermore, they 
are corrected for Galactic extinction taking the MW extinction 
law of Cardelli et al. d 19891 1 and the value of Galactic reddening 
E(B - V) from Schlegel et al. (I1998I I. Finally, they are mapped 
to the rest-frame, and smoothed to a similar resolution as for the 
FDF spectra (rest-frame 5 to 8 A). 

The CDFS was imaged with the HST ACS camera in four 
different filters as part of the ACS GOODS legacy programme 
(Giavalisco et al. 2004). We make use of the F775W band im- 
ages, where a lOcr limit of 26.5 mag was reached. For the data 
reduction, we refer the reader to Giavalisco et al. (i2004j . 

2.3. The GDDS subsample 

From the / and /Ts-limited Gemini Deep Deep Survey (GDDS, 
Abraham et al. 120041 ). we consider eight additional star-forming 
galaxies at 1.5 < z < 2.3 all of them having R < 24.7 and A's < 
21.0. The limited redshift range is mainly due to the relatively 
narrow wavelength range (typically from 5500 to 9200 A) of the 
red GMOS grating used in the GDDS. The GDDS subsample 
makes half of the total sample for 1.5 < z < 2.3 and consists of 
galaxies with a prominent, young stellar population (class '100', 
see Abraham et al. 120041 ). 

The basic spectra of Abraham et al. (2004') have been re- 
duced in the same way as described in Sect. l2.2l in order to meet 
the same standard as the spectra of the FDF and K20 subsamples 
(taken with VLT FORS). 

2.4. The reference sample of nearby starbursts 

We make use of the same reference sample of 24 local (z < 0.02) 
starburst galaxies selected by 'NP05' For these starbursts, low- 
resolution lUE spectroscopy is available for the wavelength 
range between 1 150 and 3350 A. 

Sample selection and data reduction are described in lNP05l 



3. Data analysis 

The shape of the extinction curve at rest-frame UV wavelengths 
is constrained from a suitable parametric description of the UV 
SED of a galaxy. This method was introduced by NP05 ' and 
is briefly described in Sect. 13.11 Fundamental properties of a 
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galaxy, like star-formation rate (SFR) and stellar mass, are es- 
timated by model fitting the optical spectra (Sect. [J!2l and |3.3t . 
Morphological information comes from model fits of the 2- 
D surface brightness distribution of individual galaxies with 
diffraction-Hmited HST- ACS imaging (Sect.[T4l). 

3.1. A parametric description of the UV continuum 

INP05 [ introduced a parametrisation of the rest-frame UV SED of 
a galaxy based on five power-law fits to different sub-regions of 
the UV continuum of the same form as in Calzetti et al. (.19941 ). 
Narrow wavelength regions affected by the presence of strong 
spectral lines are excluded from the fitting procedure for all pa- 
rameters (cf. Calzetti et al. 1 19941 Leitherer et al. 2002b). 

Two parameters are particularly suitable to constrain the 
slope of the extinction curve in the far-UV and the pres- 
ence/absence of the broad absorption feature centred at 2175 A. 
The first parameter characterises the apparent strength of the UV 
bump and is called 734. It is the difference between the contin- 
uum slopes measured at 1900 - 2175 A (73) and 2175 - 2500 A 
(74), respectively. A value of 734 ~ 1 indicates the absence of the 
2175 A feature. This is consistent with the presence of SMC-like 
dust in the ISM of a galaxy. Conversely, 734 < -2 points to an 
extinction law which exhibits a significant UV bump. The sec- 
ond parameter gives the amount of reddening in the UV, accord- 
ing to the SHELL models of radiative transfer by Witt & Gordon 
( 2000) . In fact, the relation between amount of dust or opacity 
and reddening is model dependent (e.g., Witt & Gordon 120001 ). 
These SHELL models with SMC dust are suitable to describe 
dust attenuation in local and high-redshift starbursts (Gordon et 
al. 119971 Vijh et al. 120031 ). Here the parameter estimating red- 
dening in the UV is defined as the continuum slope measured at 
1750 - 2600 A, with the exclusion of the range 1950 - 2400 A. 
This proxy for the UV reddening by dust is called /3b since it re- 
places the measure of UV-reddening /3 adopted by NP05 In fact, 
/? is determined at 1250 - 1750 A and, thus, is not measurable in 
optical spectra of objects at z < 2. 

In order to establish /3b as a suitable alternative to /?, we 
compare the values of these two parameters determined from 
synthetic SEDs (Fig. 13). We make use of different models com- 
bining stellar population evolutionary synthesis (Maraston l2005l ) 
and dust attenuation. For data-consistent ages of the stellar pop- 
ulations and metallicities (see Sect. |4.4t , the variations in /3b 
are only of the order of 0.1. Therefore, in Fig. [3] we assume 
the same properties for the stellar populations, i.e.: a continu- 
ous star-formation over lOOMyr at a constant rate, a standard 
Salpeter ( 1955 ) initial mass function (IMF), and solar metallicity 
Zq. Conversely, dust attenuation follows different prescriptions 
(see figure caption). In general, all models show that /3b and /3 in- 
crease together when the opacity increases, with the only excep- 
tion of models implementing the Witt & Gordon (120001 ) SHELL 
configuration with a two-phase, clumpy medium and MW-type 
dust. In this case, /3b stays constant (~ -2.6) whereas f) increases 
when the opacity increases. This behaviour is due to the fact that 
in the spectral regions where /3b is defined, attenuation increases 
at the same rate with increasing opacity. The mapping of /? into 
/3b is model dependent, of course. An overall monotonic relation 
between Pb and /3 holds even when more complex models assum- 
ing an age-dependent extinction are considered (see Sect. 14.3b . 
The locus of these models in the /3-/3b plane only moderately de- 
pends on the dust distribution, the dust-clearing time scale, and 
the star-formation history. For instance, for composite Witt & 
Gordon (i2000j models with LMC-type dust and age-dependent 
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extinction we typically find A/Sb to be less than 0.2 for constant 
p. In conclusion, we can take Pb as an alternative proxy for the 
UV-continuum reddening, though y6b has a lower dynamic range 
than p. 

Furthermore, the choice of simple but not simplistic radiative 
transfer models to illustrate the relation between ySb and [} is sup- 
ported by the fact that most of the data points occupy the region 
delimited by these models in Fig. Simultaneous measure- 
ments of y6b and (3 were possible for 38 galaxies at 1 .8 < z < 2.5 
out of the FDF and K20 subsamples, and for the 24 local star- 
burst galaxies. Part of the scatter of these data points reflects the 
dependence of ySb and /? from the dust properties of individual 
galaxies, as described by the models. Measurement errors in ySb 
contribute to this scatter as well. They are caused by the presence 
of residuals of night sky lines in the optical spectra obtained for 
the high-redshift galaxies, or by the uncertainties in the correc- 
tion for Galactic extinction and the combination of the spectra 
of the lUE short- and long-wavelength channels for the nearby 
starbursts (see lNP05l l. 

3.2. Star-formation rates 

For galaxies at z > 1.5, the available optical spectra map rest- 
frame UV wavelengths only. These spectra are contrasted to a 
suite of synthetic SEDs across the mapped wavelength domain 
excluding the range 1950 - 2400 A and intrinsic (i.e., corrected 
for attenuation by internal dust) SFRs are determined as a best-fit 
solution. In particular, we build a grid of Maraston (2005 ) mod- 
els with constant SFR0, an age between 10 and lOOOMyr, and 
metallicity equal to 0.5, 1, or 2 times solar (Z©). The ensuing 
SEDs are attenuated according to the Calzetti et al. (i2000i ) law, 
where the optical colour excess E{B-V) is the only free parame- 
ter. This choice is dictated by consistency with the literature and 
simplicity. Hence, the best-fit values of the intrinsic SFRs have 
to be taken "cum grano salis". In fact, on one hand, uncertainty is 
contributed by the unknown, true attenuation function that char- 
acterises individual objects, as discussed in Sect.|4] On the other 
hand, the data only roughly constrain age and metallicity of the 
models. Thus, the degeneracy between the analogous effects of 
age and dust reddening on the UV SED is not broken. 

For galaxies at z < 1.5, the available optical spectra do not 
map wavelengths ^ 1500 A. Therefore, we extrapolate the rest- 
frame UV SEDs of these galaxies down to the Lya line (1216 A) 
by fitting models in the UV domain excluding the range 1950 
- 2400 A. To this purpose, it is sufficient to adopt a Maraston 
( I2005I I model with age of lOOMyr, solar metallicity, standard 
Salpeter IMF and constant SFR. The ensuing SED is attenu- 
ated according to the Calzetti et al. (2000) law as before. The 
goodness of this procedure is evaluated by comparing the UV- 
continuum slope parameters for fitted synthetic rest-frame UV 
spectra and available observed spectra mapping the rest-frame 
UV down to Lya (see Sect. |2|. As a result, these fitted syn- 
thetic spectra exhibit bluer UV continua than the real spectra 
for /Imin > 1500 A. Hence, the original fitted values of E(B - V) 



^ We note that the distribution of nearby starbursts in Fig. [3] shows 
several cases of significant displacement from the locus expected from 
the Calzetti law (Calzetti et al. 1994 2000). This must not be surprising 
since the Calzetti law is a polynomial fit to the data, and, thus, expresses 
an average behaviour of dust attenuation in nearby starbursts across the 
whole wavelength range of applicability. 

' Models with a SFR that declines exponentially as a function of 
time were also explored but they were poorly constrained by the data. 
Therefore they will not be discussed hereafter. 
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Fig. 3. Comparison of the two UV-continuum slope parameters 
P (1250 < A< 1750) and Pt, (1750 < A< 2600) obtained from 
models combining stellar population evolutionary synthesis and 
dust effects as a function of the opacity of the dusty ISM. Dust 
attenuation is described either as for the empirical Calzetti et al. 
(I2OOOI 1 law (long dashed lines and circles) or as for the Witt & 
Gordon (2000) radiative transfer models for the SHELL configu- 
ration of dust and stars and a two-phase, clumpy ISM with SMC- 
(solid lines), LMC- (long dashed lines), or MW-type dust (dash- 
dotted lines). Furthermore, we make use of a screen model plus 
the extinction laws of the SMC bar (dotted lines and lozenges), 
the LMC 2 supershell (dashed lines and lozenges), the LMC av- 
erage (dashed lines and squares), and the MW average (dotted 
lines and squares) (Cardelli et al. 119891 Gordon et al. 2003). 
Symbols are plotted in intervals of AE(B- V) = 0.1. In all cases 
we adopt a Maraston ( I2005I I model with a standard Salpeter IMF, 
constant SFR, age of lOOMyr, and solar metallicity. The dia- 
gram also shows parameter values measured for 24 local star- 
burst galaxies (crosses), 33 FDF galaxies at 1.8 < z < 2.5 (filled 
lozenges), and five K20 objects at 1.9 < z < 2.3 (filled circles). 



in the Calzetti law were corrected as a function of the measured 
value of Pb ({AE{B - V)) - 0.04) and new extrapolations were 
obtained. We estimate that the uncertainty in the narrow-band, 
rest-frame UV luminosity at 1500 + 20 A, L1500 (see Noll et al. 
I2004I I, associated with this extrapolation is about 0.04 dex. 

As a consistency check, we estimate SFRs from the [Oil] 
emission doublet at 3727 A (e.g., Kennicutt |1998l Kewley et al. 
2004) for galaxies at z < 1 .5. We adopt the calibration of Kewley 
et al. (2004, see their eq. 19), assuming a nearly solar metallicity. 
We make the assumption that the values of the "instantaneous" 
SFR obtained from the [O II] line emission are statistically the 
same as the values of the SFR averaged over lOOMyr obtained 
from the rest-frame UV continuum. Fig.|4]shows that these two 
sets of estimated SFRs agree, if the stellar continuum-to-nebular 
line-emission reddening ratio is equal to 0.65 (cr = 0.08), the 
scatter in the SFR differences being of only 0.2 dex. This red- 
dening ratio is close to those estimated by Calzetti et al. (20()0]l 
and Fernandes et al. ( l2005l l. i.e., E(B-V),t^rs/E{B-V)g^^ ~ 0.44 
and 0.5, respectively. It is significantly less than the value of 1 
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Fig. 4. Comparison of SFR (UV) derived from model fits to the 
UV/optical spectra and ^/^/^([OII]) derived from [Oil] fluxes 
for the sample FDF (lozenges) and K20 galaxies (circles) at 1 < 
z < 1.5. Objects with 734 < -2 are marked by filled symbols. 
The 5F/?([OII]) were dust corrected using the Calzetti law as 
for the S FR([JW). The best agreement between these values is 
obtained for a stars-to-gas reddening ratio E{B - y)stais/£^(5 - 
^)gas = 0.65 and is reproduced here. Errors in the 5 FT? ([Oil]) 
only express uncertainties in the [O II] flux. 



often assumed in the literature (e.g., Pierini et al. 120051 Erb et al. 
l2006bD . 

3.3. Stellar masses 

Total stellar masses are taken from the Uterature. They were de- 
rived following status-of-art techniques that offer consistent re- 
sults within the uncertainties of the photometric data and model 
assumptions. They are also consistent within the uncertainties of 
a few times 0. 1 dex with the masses determined by our SFR es- 
timates (see Sect. |3.2| l. For the FDF and K20 subsamples, stellar 
masses come from Drory et al. (2005). These authors derived 
the most likely, total stellar masses and their uncertainties after 
fitting a large suite of two-component models to the FDF (Heidt 
et al. I2OO3I 1 and CDFS multi-wavelength photometry (Salvato 
et al., in prep.). Each model consisted of a combination of two 
composite stellar populations: one given by a young (100 Myr) 
burst with constant SFR, the other by an older (> 500 Myr) burst 
with an exponentially declining SFR. Both components of the 
stellar continuum emission were attenuated independently ac- 
cording to Calzetti et al. (|2000l l. For the GDDS subsample, total 
stellar masses come from Glazebrook et al. (2004). These au- 
thors derived the most likely, total stellar masses and their un- 
certainties by evaluating synthetic SEDs against the GDDS pho- 
tometry of a galaxy to determine the mass-to-K-band luminosity 
ratio, and hence the mass. Two-component models plus a stan- 
dard recipe for attenuation were adopted in a way analogous to 
Fontana et al. (|2003] ) and Brinchmann & EUis (i2000i) . 



3.4. Morphological parameters 

The package GIM2D (Simard et al. [T999t is u sed to fit point- 
spread function (PSF) convolved Sersic ('1968') profiles to the 
two-dimensional surface brightness of each object, down to a 
Umit of F814W = 24 (FDF) and FllSW = 24.5 (CDFS). 
The Sersic model contains seven free parameters: total intensity, 
scale-length of the semi-major axis, ellipticity e, position angle, 
index nser, and model x,y centring. The most valuable parame- 
ter is Hsei, which determines the shape of the profile. For each fit, 
GIM2D automatically determines the initial values and limits for 
the parameter space to be explored using a segmentation mask 
for deblending multiple objects. This mask is provided by the 
standard pipeline for the automated source extraction and pho- 
tometry SExtractor (Bertin & Arnouts 1996). The PSFs used to 
convolve the model profiles were obtained for each tile by stack- 
ing about 10 high S/N isolated stars. We put extensive eff'orts in 
understanding and quantifying how various kinds of systematics 
(e.g., PSF variations over the field of view or dependences on 
the S/N ratio) could affect the modelling results. To this purpose, 
simulated images with the same characteristics of the true ones 
were produced and analysed in the same way as for the ACS 
images. GIM2D also computes the two indices /?t and Rp, from 
the thumbnail residual image. These indices provide an estimate 
of the overall smoothness of a galaxy image with respect to the 
fitting model. In other words, they estimate the residual substruc- 
ture like spiral arms in nearby late-type galaxies (see the seminal 
study of Elmegreen et al. 1992), peculiarity/asymmetry in the 
distribution of giant star-forming regions, or interaction/mergers 
in high-redshift galaxies (e.g., Schade et al. 1995) . 

Furthermore, we characterise the (potentially complex) mor- 
phologies of our galaxies by means of the model-independent 
parameters concentration C and asymmetry A. The C-A method 
was developed in the mid-nineties by Abraham et al. (1994] 
[1996). Subsequent work s (e.g Wu ,1999. Consehce et al. .20001 
2003; Menanteau et al. I2006I I have shown that a better mor- 
phological classification is obtained by choosing an image pivot 
point which minimises the measured asymmetry. We make use 
of the CAS parametrisation as proposed and described in detail 
by Conselice et al. ( 120001 |2003l l. Here it is important to say that 
early-type galaxies have larger concentration and lower asym- 
metry indices than later ones. 

4. Results 

Hereafter we show that dust with a UV bump can be present 
in massive, star-forming galaxies at 1 < z < 2.5 (Sect. 14. 11 1. 
Sect. |4.2| contains reasons why sample selection criteria can hin- 
der the inference of the presence of the dust components respon- 
sible for the absorption excess at 2175 A (i.e., the carriers of 
the UV bump) from a similar spectral analysis of high-redshift 
galaxies. We interpret our results in Sect. 14.31 Then we inves- 
tigate how the presence of different extinction curves in high- 
redshift galaxies depends on global galaxy properties such as 
SFR and total stellar mass (Sect. I4.4l i. Finally, we discuss the 
effects of topology and metal enrichment of the ISM on the ob- 
served strength of the UV bump (Sect. |43] l. 

4. 1 . Dust with a pronounced UV bump already exists at high 
redshift 

Fig. |5] shows the proxy for the strength of the 2175 A feature, 
734, versus the proxy for reddening in the UV, fib, for three dif- 
ferent redshift bins, i.e.: 1 < z < 1.5 (40 objects), 1.5 < z < 2 
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Fig. 5. Plot of the proxies for the 2175 A feature (734) and red- 
dening in the UV (j6/,) for the FDF (lozenges) and K20 (circles) 
galaxies listed in Table [1] for which both parameters were deter- 
mined. The panels show the distribution of objects at 1 < z < 1 .5 
(top), 1.5 < z < 2 (middle), and 2 < z < 2.5 (bottom), to- 
gether with that of the 24 local starburst galaxies (crosses) as a 
reference. All panels reproduce model tracks obtained for dif- 
ferent descriptions of dust attenuation as a function of E(B - V) 
(see Fig.[3j. We recall that more negative values of 734 indicate 
stronger 2175 A absorption features. 
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Fig, 6. Comparison of /Tj -weighted composite spectra of 1 < z < 
1.5 galaxies with R < 24 exhibiting 734 > -2 (blue or dark 
grey) and 734 < -2 (red or grey), respectively (top). The ratio 
of both composites, normalised at 2400 - 2570 A, is also shown 
(bottom). 



(10 objects), and 2 < z < 2.5 (38 objects). Tracks for the models 
described in Sect. 13. H are also plotted. The comparison between 
the distribution of the observed quantities and the expected ones 
clearly shows that dust with a UV bump and, thus, different ex- 
tinction curves (see Sect. 14.3b do exist at high redshift. Now this 
holds also for 1 < z < 1.5 (cf.lNPOSl). 

Galaxies at 1 < z < 1 .5 exhibit the widest range in the ob- 
served strength of the UV bump (-8 < 734 < 1). According 
to the models described in Sect. 13. II dust with a significant UV 
bump (i.e., 734 < -2) characterises 58% of the 19 /Tj-selected 
K20 galaxies and 19% of the 21 /-selected FDF galaxies at 
these redshiftfl This drop is partly caused by selection effects 
(see Sect. 14. 2t . In fact about 29% of the galaxies in the original 
FDF photometric sample with 1 < Zphot < 1 5 and R < 24 have 
K, < 20 (Heidt et al. l2U03l Gabasch et al. l2004] l. The fraction of 
FDF galaxies with 734 < -2 reaches 33%, if only the six objects 
with < 20 are considered. This suggests that dust with a sig- 
nificant UV bump can be present in at least half of the /Tj -bright, 
star-forming galaxies at 1 < z < 1.5, whatever the selection cri- 
terion. 

Furthermore, Fig.|5]shows that the 2175 A absorption feature 
tends to be associated with larger amounts of reddening at UV 
wavelengths. This confirms the result of |NP05l Interestingly the 
largest amounts of reddening appear in the lowest redshift bin. 
When we combine the K20 and FDF objects with 1 < z < 1.5 
to determine the A's-weig htecfl mean value of 734, we find that 
(734) = -0.44 + 0.18 if ySb < -1.5, whereas (734) = -3.48+0.51 
if Pb > -0.5. Consistently, in the same redshift bin, about 80% 
of the /Ts-bright galaxies with /3b > -0.5 exhibit evidence of dust 



"* We note that a MW-like extinction law is directly suggested only 
forCDFS-00525. 

' For subsamples the weight of A's-bright objects differs from the 
value of 29% for the whole sample owing to the different distributions 
of the A^s-bright and A's-faint galaxies as a function of individual selec- 
tion quantities. 
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Fig. 7. The distribution of 734 values as a function of redshift 
for the full sample of 108 massive, star-forming galaxies at 
1 < z < 2.5 selected from the FDF Spectroscopic Survey 
(lozenges and x), the K20 Survey (circles and +), and the GDDS 
(asterisks). Crosses and asterisks mark galaxies with no possible 
determination of the UV continuum slope /3h owing to a limited 
spectral coverage. Filled symbols indicate objects with 734 < -2, 
i.e., with an extinction curve exhibiting a significant UV bump. 



with a UV bump, but only about 40% of those with /3h < -0.5 
do so. 

The evidence of a UV bump (734 < -2) in Fig.|5]is made 
clearer by the comparison between the /Tj-weighted, composite 
spectra of galaxies at 1 < z < 1.5 with R < 24 and 734 > -2 or 
734 < -2 (Fig.|6]). An analogous proof was given by ,NP05T or 
their sample of FDF UV-luminous galaxies at 2 < z < 2.5 (also 
reproduced in Fig.|5]). Owing to the different mapping of the rest- 
frame UV domain, the composite spectra in Fig.|6]show remark- 
able differences in the strength of the Balmer/4000 A break and 
in the Mg II doublet at 2800 A. The origin of these differences is 
discussed in Sect. 14.41 



4.2. Selection effects on the detection of ttie UV bump 

In Fig.|5] the few galaxies at 1.5 < z < 2 do not seem to possess 
significant 2175 A features in their spectra, which is at odds with 
the UV-bump detections at lower and higher redshifts. Hence, 
these objects may not represent the full population of massive, 
star-forming galaxies at 1.5 < z < 2, though their statistics is 
poorer than for objects in the other two redshift bins. Fig. [8] 
shows that the sample galaxies at 1.5 < z < 2.2 (all exhibit- 
ing 734 > -2) have rest-frame luminosities at 1500 A that are 
on average three times higher than those of the galaxies with 
734 < -2 at 1 .0 < z < 1 .5 but slightly less luminous than those 
with an evident UV bump at 2.3 < z < 2.5 (34.02 + 0.06 versus 
34.16 + 0.04 in log Lisoo - in [W/A], respectively). This sug- 
gests that our selection can miss galaxies with a low, intrinsic 
UV-luminosity and a significant UV bump, such as those typi- 
cally found at 1.0 < z < 1.5. However, it could also miss objects 
with a high, intrinsic UV-luminosity and 734 < -2, such as those 
found at 2.0 < z < 2.5. This can be understood if the presence of 
a significant UV bump is indeed associated with a larger amount 
of attenuation at UV wavelengths, as suggested by the model in- 
terpretation of Fig.|5] In fact, the sample galaxies at 1 .5 < z < 2.2 
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Fig. 8. Luminosity at 1500 A versus redshift for 88 galaxies of 
the FDF Spectroscopic Survey (lozenges) and the K20 survey 
(circles). Galaxies with 734 < -2 are marked by filled symbols. 



are mostly selected in / band and exhibit R < 24 (cf. Sect. |2]l. 
Now, we note that the R broad-band filter maps exactly the spec- 
tral region across the UV bump for 1.5 < z < 2.2. As a check of 
this hypothesis, we plot the distribution of 734 as a function of 
redshift for the full sample of 108 massive, star-forming galaxies 
at 1 < z < 2.5. With respect to Fig.|5] Fig. [T] contains additional 
20 galaxies with no measurement of /3\,. Among the 26 galaxies 
at 1 .5 < z < 2.2, there are only two with 734 < -2 against an ex- 
pected number of 8+3 (corresponding to a fraction of about 30% 
as from Sect. [4311. For one (CDFS-00295) the determination of 
its spectroscopic redshift is uncertain. The other (SA15-4762) 
has R - 24.4, i.e., it is just fainter than the limiting magnitude 
(R = 24) of the FDF and K20 subsamples. We conclude that se- 
lection effects can remove from the resulting sample very dusty 
galaxies with dust similar to the LMC or MW at 1.5 < z < 2.2. 

4.3. Extinction curve and dust/stars configuration 



In the previous analysis as well as in'NP05' it was assumed that 
dust attenuation in massive, UV-luminous galaxies at high red- 
shifts is described by the Witt & Gordon (2000) Monte Carlo 
calculations of radiative transfer for the SHELL configuration. 
In particular, a two-phase, clumpy, dusty medium is assumed 
to surround all stars, whatever the ages of the stellar popula- 
tions. If this assumption holds, the interpretation of our results 
is straightforward: different mixtures of dust produce extinction 
curves that are intermediate between the SMC and LMC average 
extinction curves. 

However, the previous assumption on the dust/stars configu- 
ration can sound too simple and, thus, be challenged. Therefore, 
hereafter we discuss two dust/stars configurations that differ sub- 
stantially from each other and from the SHELL configuration in 
Witt & Gordon ( 120001 ). In the first case, we consider a disc ge- 
ometry together with a different distribution of stars with respect 
to the one of the dust according to the different ages of the stellar 
populations. This effect is known as the "age-dependent extinc- 
tion" (Silva et al. 11998 ; Tuffs et al. 2004; Pierini et al. |2004t 
Panuzzo et al. | 2007| l. In the second case, we assume that any 
of our sample galaxies can be ideally pictured as a system where 
stellar populations younger than a certain age (i.e., the dust clear- 
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Fig. 9. Comparison among models where stars and dust are dis- 
tributed in a disc. First, dust attenuation is described as in Pierini 
et al. (I2004I I. assuming inclinations of (solid line), 70 (dashed 
line), and 90 deg (dotted line) and continuous star-formation 
phases of 600 (small open circles) or 2000 Myr (small filled cir- 
cles). The V-band opacities plotted are 0, 0.5, 1, 2, 4, 8, and 16. 
Second, we show Sc-galaxy GRASIL models (Silva et al. 1 19981 1 
for inclinations of 0, 75, and 90 deg (crosses, from the left to the 
right, respectively). 



ing time t^) are enshrouded by dust, while older ones live in a 
dust-free environment. 

We consider two different sets of models describing dust at- 
tenuation for a disc geometry, namely the "Sc model" from the 
GRASIL library (Silva et al. |1998| ) and the d isc mo dels with a 
two-phase, clumpy, dusty ISM in Pierini et al. (12004 1) . In general, 
the GRASIL radiative transfer models describe the fact that new 
stars are born inside molecular clouds and at later times either 
get rid of their parent environment or move out of it. The radia- 
tive transfer models described in Pierini et al. (2004) account for 
the fact that younger stellar populations are more attenuated by 
dust than older ones by setting the scale-height of the distribution 
of photons equal to lower values the shorter the wavelength. In 
both cases a MW-like extinction curve is assumed. Three values 
of the inclination of the disc are considered in each set of mod- 
els, namely (i.e., a face-on view), 75 (or 70), and 90 (i.e., an 
edge-on view) degrees. Furthermore, six different values of the 
opacity are considered in the second set of models. In both cases, 
intrinsic SEDs are computed for a constant SFR, solar metallic- 
ity, and ages of 0.6 or 2Gyr. As Fig.|9]shows, the models do not 
span a domain as wide as that of the data in the fib-ju plane 
(cf. Fig. Is). In particular, they are confined to -2.5 < 734 < 1 
although they assume an extinction curve similar to the average 
MW one. Differences do exist however The GRASIL Sc models 
span a limited range in fih as a function of the inclination of the 
disc with respect to the disc models described in Pierini et al. 
( I2004I I. once the opacity is fixed. Assuming the presence of dust 
with a UV bump and a smaller scattering albedo for a shorter 
wavelength (except for the UV-bump range) would broaden the 
range in fi\, spanned by the previous models but would not move 
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Fig. 10. Comparison between the distribution in the fih-jM plane 
of observed galaxies at 1 < z < 1.5 and models where the star- 
formation activity takes place in "macro-regions" until dust is 
cleared out there. Dust attenuation is described as in the Witt 
& Gordon (2000) SHELL configuration, where the ISM is ei- 
ther two-phase, clumpy (thin lines) or homogeneous (thick lines) 
and contains average LMC-type dust. As for (see Sect. 14. 31 ). it 
is set equal to either 100 Myr (dash-dotted and solid lines) or 
400 Myr (dashed lines). Star-formation goes on for a time equal 
to the dust clearing time fc (solid, black lines; cf. Fig. |5]) or to 
1 Gyr. For the latter the total stellar mass produced by the extin- 
guished stellar population Mext is fixed to be half (green or light 
grey lines) or twice (red or dark grey lines) the mass produced 
by the unextinguished population Munext- These constraints lead 
to S FRsxi/S FR^nfijii between 0.75 and 18. Three of four star- 
formation scenarios prescribe the instantaneous SFR at t{z) to be 
above the SFR averaged over the age of the system. 



any of them to the region of the highest values of y6b and the low- 
est values of 734 (cf. fig. 1 1 in Inoue et al. 2006). In conclusion, 
the comparison of observed values and predictions from a more 
physical dust/stars configuration for a disc geometry leads us to 
the conclusion that especially those galaxies with 734 < -2 seem 
to have dust with a significant UV bump that is partly distributed 
out of the disc mid-plane or at least embeds the stellar popula- 
tions responsible for the UV emission. This holds whatever the 
behaviour of the dust scattering albedo. 

Alternatively, we assume that in our high-redshift galaxies 
new stars are born in "macro-regions" where the dusty ISM is 
cleared out after a time t^, which ends the star-formation activ- 
ity there. For illustrative purposes only, we make the simplistic 
case that the stellar populations younger than the dust clearing 
time are screened by dust while the older ones do not suffer from 
extinction at all. As for the value of t^, we do not take the time 
scale of the dispersal of the parent molecular cloud by an in- 
dividual supernova, as in Silva et al. ( 119981 ) or in Tuffs et al. 
(2004), which is of the order of a few Myr (e.g., Leitherer et 
al. 2002a). Conversely, we picture dust clearing as a collective 
phenomenon and set equal to 100 or 400 Myr Dust attenu- 
ation is described as for the SHELL configuration in Witt & 
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Gordon ( 120001 1. assuming either a homogeneous or two-phase, 
clumpy ISM with LMC-type dust. Only values of opacity giv- 
ing a maximum attenuation at 1500 A of 5 mag are considered. 
Furthermore, we consider Maraston ( |2005l l stellar population 
synthesis models with maximum ages of 1 Gyr and fixed, so- 
lar metallicity (see Sect. 14.4. It . The total stellar mass produced 
by the extinguished stellar population Mgxt is fixed to be half 
or twice the mass produced by the unextinguished population 
MuiKxt- The SFR is, therefore, changed (divided by a factor from 
0.75 to 18) after a time equal to t^. Fig.llOlshows a comparison 
of the previous models and the data in the ySb-7'34 plane, only for 
galaxies at 1 < z < 1.5. Remarkably, this set of models spans 
the region defined by the observed measurements. For low opac- 
ities the light from extinguished, young stars dominates the UV 
spectra and the observed strength of the UV bump increases with 
UV reddening. For higher opacities the contribution of the un- 
obscured, older stars becomes more and more important till 734 
stalls or the strength of the 2175 A feature even decrease^ The 
latter trend is more striking for galaxies with short dust-clearing 
time scales. The most negative 734 in our data can only be mod- 
elled assuming little impact of nearly-unobscured intermediate 
age stars on the UV SEDs and the presence of an almost perfect 
dust screen for the obscured stellar populations. LMC-type dust 
or dust with even stronger UV bumps are required. 

In conclusion, age-dependent extinction scenarios increase 
the probability of the presence of dust with a significant UV 
bump in our sample galaxies. 

4.4. Extinction curve and global galaxy properties 

The carriers of the UV bump exist in galaxies belonging to dif- 
ferent subsamples. 2175 A features traced by 734 < -2 are found 
for about 38% of the FDF galaxies at 2.3 < z < 2.5. They exist 
as well in the spectra of about 25% of the A's -weighted combined 
subsample of FDF and K20 galaxies with /J < 24 at 1 < z < 1.5. 
This fraction jumps to about 52% if one considers only TSTs-bright 
galaxies at these redshifts (Sect. 1431 ). On the other hand. Fig. [8] 
shows that galaxies at 1 < z < 1.5 exhibit lower luminosities at 
1500 A than those at 2 < z < 2.5, in general. For objects with 
evidence of dust with a UV bump, L1500 drops by a factor of 4 
- 5. This suggests that the reasons for the presence and/or sur- 
vival of the carriers of the UV bump in high-redshift galaxies are 
different, at least in part, at different redshifts. 

Hereafter we investigate if, more generally, the extinction 
curve depends on global characteristics of a high-redshift galaxy 
like properties of the stellar populations, SFR, stellar mass, and 
morphology. We will consider six subsamples, which are clas- 
sified by redshift, UV-continuum slope, and observed strength 
of the UV bump. In particular, we will compare a sample at 
1 < z < 1.5 (/Ts -weighted, see Sect. 14.11 ) and one at 2 < z < 2.5 
(only FDF galaxies). Both samples are divided into three sub- 
samples each. First, we define "blue" galaxies, i.e., galaxies with 
blue UV continua. For 2 < z < 2.5, we take the same classifica- 
tion as introduced by'NP05', i.e., the "blue" galaxies (16 objects) 
have /3 < -0.4. For 1 < z < 1.5 we take /3b < -1.5 (13 ob- 
jects). The remaining, so-called "red" galaxies are divided into 
subsamples for 734 > -2 (14 objects at 2 < z < 2.5 and 9 ones 
at 1 < z < 1.5) and 734 < -2 (13 and 9 objects, respectively), 
since galaxies with relatively strong 2175 A features are mostly 
"red", i.e., /3 > -0.4 and /3h > -1.5, respectively. We choose 
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* An analogous decrease of the observed strength of the UV bump 
(when the fraction of less obscured stars older than but still emitting 
UV photons increases) was already found by Panuzzo et al. ( |2007| l. 



Fig. 11. B aimer-break proxy C(36-39) versus reddening measure 
ySb for the sample FDF (lozenges) and K20 galaxies (circles) with 

1 < z < 1.5. Objects with 734 < -2 are marked by filled symbols. 
Different Maraston ( 2005) models for continuous star-formation 
and solar metallicity are shown for estimating the age of the stel- 
lar population. The dashed lines indicate the effect of continuum 
reddening by a Calzetti et al. (|2000|) law on ySb and C(36-39) for 
ages of 10, 100, lOOOMyr. The small filled circles are plotted in 
intervals of AE{B - V) = 0.1. Attenuation via the Calzetti law 
produces a stronger dependence of C(36-39) on /3h than most 
other attenuation laws used in this study (see Fig. [3]). Therefore, 
the estimated ages would typically increase by a factor of about 

2 using SMC and LMC extinction curves. The dust-free cases 
for the three considered ages are indicated by the dotted lines. 



y6b = -1.5 for 1 < z < 1 .5 since it is consistent with the discrim- 
inating value of y6 for the higher-redshift galaxies. In fact, for 
those objects at 2 < z < 2.5 with a simultaneous measurement 
of j6 and /3b we find /3 - /3b ~ 1 . Comfortably, an analogous dif- 
ference is obtained for the reference sample of nearby starbursts. 
By chance, ySb = -1-5 divides the FDF subsample of galaxies at 
1 < z < 1.5 into blue and red subsamples of similar sizes (11 
and 10 objects, respectively). 

4.4.1. Stronger 2175 A features in more "mature" galaxies 

Fitting the available spectra with the models described in 
Sect. 13. 2| gives ages of less than lOOMyr (for a constant SFR) 
for the majority of the sample galaxies. This simply confirms 
the star-forming nature of the sample galaxies. Ages of more 
than 1 Gyr were determined for high-redshift star-forming galax- 
ies like ours by Erb et al. ('2006b) and Glazebrook et al. ( 2004l l. 
based on broad-band SEDs mapping the full rest-frame optical 
domain. 

For galaxies at z < 1.3, the spectra map also rest-frame 
optical wavelengths, reaching the region of the Balmer/4000A 
break, where the index D4000 is defined (Hamilton 119851 
Balogh et al. 1999; Kauffmann et al. 2003). This spectral in- 
dex is particularly useful to constrain ages > 1 Gyr, also at high 
redshifts (e.g., Le Borgne et al. |2006I ). We adopt the definition 
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of D4000 by Balogh et al. d 19991 1. where two fluxes (per unit 
of frequency) are measured in the wavelength windows 3850 - 
3950 A and 4000 - 4100 A, bracketing the Balmer/4000 A break. 
We measure a typical flux ratio of about 1.1, which translates 
into an average age of 50^^^**Myr for Maraston (2005 ) mod- 
els with constant SFR and attenuated according to the Calzetti 
law, which provides lower age limits (see caption of Fig. fTTT l. 
In the complementary dust-free case, this estimate increases to 
250^j^Q Myr. For more complex star-formation histories, includ- 
ing both a monotonic behaviour of the SFR as a function of time 
and a random burst of star-formation (cf. Kauffmann et al. 2003 ), 
a value of D4000 - 1.1 corresponds to a median age of the mod- 
els equal to 630Myr (A. Gallazzi, priv. comm.). The ages of the 
latter models have a large scatter though, ranging from 126 Myr 
to 4 Gyr True characteristic ages will lie somewhere in between 
the dusty and dust-free cases. 

Characteristic ages can be determined for all objects at z < 
1.5, if one considers the Balmer-break sensitive pseudo-colour 
C(36-39), defined as the flux ratio (in magnitudes) between the 
fluxes (per unit of frequency) at 3500 - 3700 A and 3750 - 
3950 A (cf. Mignoli et al. I2005I I. In this case, we derive a mean 
age of 25^5° Myr from the comparison between data and models 
(Fig. fTTT l. If only galaxies at z < 1.3 are considered, the mean 
age increases to SO^jj] Myr This is consistent with the indica- 
tions from D4000. As also suggested by Fig. |6] the mean ages 
of galaxies with 734 < -2 are greater than those of galaxies with 
734 > -2 (50+ 15 Myr and 20+5° Myr, respectively). This trend is 
strengthened in the dust-free case, where upper limits to the char- 
acteristic age of the stellar populations are equal to 380^jgQMyr 
and 80;^3Q Myr for galaxies with 734 < -2 and 734 > -2, respec- 
tively. In conclusion, the analysis based on spectroscopic indices 
reveals the presence of intermediate-age (0.2 to 1 - 2 Gyr-old) 
stellar populations in galaxies at 1 < z < 1.5. Furthermore, it 
suggests that galaxies with evidence of a UV bump at these red- 
shifts have a larger fraction of intermediate-age stars and/or are 
older than those lacking an evident UV bump. 

Standard metallicity indicators like the R23 index (see, e.g., 
Kobulnicky & Kewley 2004) are not accessible to optical tele- 
scopes for galaxies at high redshifts. So we estimate the metallic- 
ities of galaxies at 1 < z < 1 .5 from the mass-metallicity relation 
evolution model by Savaglio et al. ( 2005), which was established 
for galaxies at z < 1. We obtain 12 + log(0/H) = 8.7 + 0.2, i.e., 
a metallicity close to Zq. We caveat the reader that metallicity 
estimates based on nebular lines overestimate the stellar metal- 
licity by up to 60% (e.g., Bresolin et al. 2005). Nevertheless, a 
solar-like metallicity is supported by a curve-of-growth analy- 
sis of interstellar metal absorption lines in 1.3 < z < 2 GDDS 
galaxy spectra (Savaglio et al. |2004| l. 

For galaxies at z > 1.5, we can use the J - K colour as a 
proxy for the characteristic age of a galaxy, since J - K brackets 
the Balmer/4000 A break for 2 < z < 4. The subsample of FDF 
galaxies at 2 < z < 2.5 exhibits {J - K) - 1.7 (ctj^k: - 0.4). 
Hence these galaxies do no meet the J - K > 2.3 criterion 
for distant red galaxies at z > 2 (DRGs, Franx et al. 120031 
van Dokkum et al. 2004), so that the bulk of their stellar pop- 
ulations is most probably not much older than 1 Gyfl (see also 
Pierini et al. 2005 ). A certain degree of evolution is required 
by the level of enrichment of the stellar populations reached by 
galaxies at 2 < z < 2.5 in the FDF subsample (see NP05'). In 
fact, for all FDF galaxies at these redshifts with available high- 

' This age constraint is compatible with the best-fit age estimates of 
three DRGs (0.3, 1.3, and 1.4 Gyr) of Kriek et al. (iSoEi . based on 
D4000 indices ranging from 1.2 to 1.4. 



S/N spectroscopy at rest-frame UV wavelengths, Mehlert et al. 
( 120021 |2006t determined typical metallicities of 0.5 - 1 Z.0, us- 
ing the stellar-wind feature CYV A 1550 and the photometric in- 
dices '1370' and '1425' (Leithereret al. |2001| l, respectively. This 
range is in agreement with the mass-metallicity relation for z ~ 2 
ofErbet al. (I2006al l. 

4.4.2. No trend between extinction curve and star-formation 
rate or stellar mass 

Table |2] lists average values of the best-fit SFR obtained for in- 
dividual galaxies (Sect. 13.21 ) that are grouped in different sub- 
samples. In general, less reddened galaxies exhibit lower SFRs 
at any redshift as a consequence of a lower value of the best- 
fit amount of attenuation. In fact, blue subsamples typically ex- 
hibit E(B - V) ~ 0.2 - 0.3, while red subsamples seem to have 
E(B -V)~ 0.4 on average, whatever the redshift (Table|2|. 

We note that the SFRs of the 2 < z < 2.5 subsamples are gen- 
erally higher by a factor of 2 to 4 than those of the corresponding 
1 < z < 1.5 subsamples. The largest differences are observed for 
red galaxies associated with a significant UV bump. In fact, red 
galaxies at 2 < z < 2.5 tend to exhibit a stronger absorption 
feature centred at 2175 A the larger the attenuation, and, thus, 
the higher the SFR. Consistently, these galaxies are much redder 
than typical Lyman-break galaxies (cf Shapley et al. 2004; Noll 
et al. 120041 ) and have higher SFRs (about 500 Mq yr ' on av- 
erage). Their SFRs exceed those estimated for galaxies selected 
by UV colour indices (e.g., Erb et al. I2006bl ) by about one or- 
der of magnitude. Daddi et al. ( 120041 ) estimated SFRs of almost 
the same order of magnitude for the K20 galaxies at the highest 
redshifts in our sample. 

In addition. Fig. [12] shows the average total stellar mass as a 
function of 734 for different subsamples. The average total stellar 
mass is equal to 1.3(+0.3) x lO'^Mo and 2.8(+0.5) x IO'^Mq 
for galaxies at 1 < z < 1.5 and 2 < z < 2.5, respectively. It 
appears to be larger in star-forming galaxies at 1 < z < 2.5 
that suffer from a larger amount of UV-continuum reddening. 
However, there is no trend between observed strength of the UV 
bump and total stellar mass, as already found bv lNP05l 

Finally, we determine specific SFRs, i.e., SFRs per unit of 
stellar mass. As shown in Table|2l the galaxies in our sample ex- 
hibit a narrow range in specific SFR {(p ~ 10' Gyr"') whatever 
the redshift. This is a rather large value, especially when consid- 
ering that (p has a range of four orders of magnitude for much 
larger samples that probe the mixture of galaxies at high redshift 
(e.g., Reddy et al. 2006 ). A plausible reason for the narrow range 
of high (p of our galaxies is that they are all close to the peak of 
their star-formation history. 

4.4.3. Extinction curve and morphology 

Table |2] lists average values of the effective radius 7?e, Sersic 
index «ser, concentration C, asymmetry A, and "dumpiness" 
Rj + Ra of galaxies grouped in the six subsamples that were 
described previously. We stress that these morphological param- 
eters were determined in different rest-frame wavelength do- 
mains from the mid-UV to the U band depending on redshift (see 
Sect. 13.4b . In general, the rest-frame UV/f/-band morphology of 
a galaxy does not appear to be directly related with the shape 
of the extinction curve at UV wavelengths whatever the redshift. 
Nevertheless, blue objects at 2 < z < 2.5 and 1 < z < 1.5 tend to 
have smaller effective radii than red objects at the same redshift. 
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Fig. 12. SFR (top) and total stellar mass (bottom) versus the 
proxy for the 2175 A feature 734 for the six subsamples defined 
in Sect. 14.41 Lozenges and circles indicate the FDF 2 < z < 2.5 
and the FDF+K20 1 < z < 1.5 galaxies, respectively. Open 
symbols mark the subsamples with low reddening (J3 < -0.4; 
/?b < -1-5), while the highly reddened galaxies are represented 
by filled symbols. Mean errors are indicated. 



whether or not a UV bump is detected in their spectra (see also 
INP05]) . 

As for the Sersic index, there is a hint that red galaxies with 
evidence of a UV bump have a lower Wser than red galaxies with 
no evidence of it in both redshift bins (see Fig.flJ]). For instance, 
at 1 < z < 1.5, the former exhibit an average Sersic index of 
0.54 + 0.25, the latter one equal to 0.98 + 0.23. Overall, the 
Sersic index seems to be larger for galaxies at higher redshifts 
(Fig.flJTl. Since the available high-resolution imaging probes the 
rest-frame UV/f/-band morphology, the Sersic profile describes 
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Fig. 13. Sersic index (top) and concentration (bottom) versus the 
proxy for the 2175 A feature 734 for the six subsamples defined 
in Sect. 14.41 Symbols are as in Fig. [12] Mean eiTors are indicated. 



the large-scale distribution of star-formation regions instead of 
the classical Hubble type. Hence, the azimuthally-averaged ra- 
dial distribution of the rest-frame UV/C/-band light traced by 
«sei is shallower for galaxies with dust producing a significant 
2175 A feature with respect to those without it, whatever the red- 
shift. Furthermore, it is more peaked at higher redshifts than at 
lower ones, despite this may be an effect of the cosmological 
brightness dimming. In fact, a faint component of a galaxy like a 
disc can fail detection more easily at higher redshifts, thus pro- 
ducing a spuriously larger value of «sei . 

In addition to the relatively large effective radii and low 
Sersic indices, the low values of concentration of the red galax- 
ies with evidence of a UV bump at 1 < z < 1.5 (see Fig. \T3[ 
suggest that most of these galaxies are large systems with shal- 
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Fig. 14. CAS parameters concentration C and asymmetry A 
(Conselice et al. 120001 12003I I for our sample of FDF (lozenges) 
and K20 galaxies (circles) at 1 < z < 1.5. Galaxies with 734 < 
-2 are marked by filled symbols. The solid and dashed lines sep- 
arate early-type, late-type, and merging galaxies by their asym- 
metry (see Conselice et al. 120031 1. 



low radial profiles at rest-frame UV/f/-band wavelengths. This is 
not surprising since most of the galaxies at these redshifts can be 
classified as late types (see Fig.[T4ll. However, the visual inspec- 
tion of the ACS images reveals that the fraction of objects with 
a shallow light profile in the rest-frame UV/fZ-band is 40 - 70% 
if /3b > -1.5, but only 20 -40% if ySb < -1.5. For red galaxies 
with a significant UV bump, this fraction becomes larger (60 - 
80%), though there are exceptions: CDFS-0271, the object with 
the strongest observed 2175 A absorption feature, appears as a 
quite compact galaxy (/?e = 2.5 kpc). 

Establishing the rest-frame UV morphology of a galaxy at 
2 < z < 2.5 is not possible from the available data, owing to the 
increased cosmological dimming. Nevertheless, we can investi- 
gate the presence of strongly distorted morphologies and/or mul- 
tiple main components. These characteristics exclude the possi- 
bility of a single object with a smooth, radial surface brightness 
profile in the rest-frame UV. While at 1 < z < 1.5 almost all 
galaxies show only one major component, at 2 < z < 2.5 45 - 
65% of the objects seem to have two or more main components. 
There is not much difference between red (fi > -0.4) and blue 
galaxies (J3 < -0.4) as for the fraction of objects with multiple 
components. However, there is a considerable difference when 
only red objects with strong 2175 A features are considered. 
For these galaxies, the fraction of objects with multiple compo- 
nents rises to 70 - 80%. This suggests that galaxies with extinc- 
tion curves exhibiting a significant UV bump at 2 < z < 2.5 
are either systems with many, large star-formation complexes 
or merging systems. This could point to the existence of an in- 
trinsic structural difference with respect to analogous galaxies at 
1 < z < 1.5, which appear as smooth, disc-like systems in the 
rest-frame UV/f/-band. 



4.5. Properties of the dusty interstellar medium 

For the subsample of galaxies at 2 < z < 2.5, the strength of the 
2175 A feature is linked to the average equivalent width (EW) of 
six prominent, interstellar low-ionisation, absorption lines Wlis 
mapped by the FDF spectra (NP05). This relation suggests that 
the presence of a significant UV bump is linked to high col- 
umn densities of (cold) gas and dust in the direction towards 
the observer, since Wus traces the covering factor of neutral-gas 
clouds (Shapley et al. 2003 ). In this scenario, the carriers of the 
UV bump are protected from strong and hard radiation fields by 
other (more robust) dust grains (e.g., Gordon et al. 2003). For 
objects at 1 < z < 1.5, there is a lack of strong, purely interstel- 
lar absorption lines at the rest-frame UV wavelengths mapped 
by optical telescopes. Nevertheless, one can investigate the most 
striking absorption features, i.e., the Fell doublet at 2600 A and 
the Mg 11 doublet at 2800 A. We find that the strength of the Fe 11 
absorption does not depend on y6b or 734, being Wpeii ~ 5 A for 
all subsamples. On the other hand, WMgii = 1.8 + 0.7 A for blue 
galaxies but rises to 5.0 + 0.5 A for red galaxies with a strong 
2175 A feature. Hence, the strength of the Mgll absorption in- 
creases together with the apparent strength of the UV bump. 
The existence of this trend is robust against the contamination 
of the measured Mg 11 equivalent width by an emission contri- 
bution from an active galactic nucleus (AGN) since the presence 
of AGN in the sample is excluded by the line profiles. 

The explanation of these results is not straightforward, given 
the multiple origin of the absorption in the Fe 11 and particularly 
the Mgll doublets. In fact, the latter feature originates in stellar 
atmospheres as well as in the ISM and, thus, depends on the frac- 
tion of intermediate-type stars (particularly A stars), the metal- 
licity and chemical composition of the stellar populations as well 
as of the ISM, and the structure of the ISM. With these caveats in 
mind, even a significant contribution of photospheric absorption 
(see, e.g., Fanelli et al. |1992| l to the strength of Mgll will not 
harm our viable interpretations, that are based on a purely inter- 
stellar origin of the Mgll and Fell lines. First of all, we note that 
Mg and Fe have similar condensation temperatures and their de- 
pletions exhibit a similar systematic trend with density (Whittet 
2003). However their relative gas-phase abundances depend on 
the gas temperature (Savage & Sembach 1996). In the warm- 
and cold-gas phases, there is less than 10% and 1% of Fe, re- 
spectively. Furthermore, Fe is depleted on dust grains about six 
times more than Mg in the cool ISM, but the Mg-to-Fe gas-phase 
abundance ratio decreases to about 2 in the warm ISM (Savage & 
Sembach 1996). Hence, a larger Mg-to-Fe equivalent width ra- 
tio can be explained independently by a greater metallicity (plus 
depletion on dust grains), as a larger fraction of the ISM being 
cool or as the ISM being particularly enriched in Mg. Sources of 
enrichment of the ISM, with a high yield of Mg but only a mod- 
est one (if any) of Fe, can be identified as the Asymptotic Giant 
Branch (AGB) stars. They dominate the overall mass-loss bud- 
get at least in the Galaxy (Whittet [20031 ). Both the presence of 
these intermediate-mass (2-5 Mq) stars and a higher metallicity 
are naturally consistent with the ages estimated for red galaxies 
with evidence of a UV bump and a larger Mg-to-Fe equivalent 
width ratio at 1 < z < 1.5 (see Sect.lSj). 



5. Discussion 

In the local Universe, well-characterised extinction curves exist 
for different environments in the Milky Way and the Magellanic 
Clouds. From the MW to the SMC, the presence of a broad ab- 



S. Noll et al.: Presence of dust with a UV bump in massive, star-forming galaxies at 1 < z < 2.5 



13 



sorption feature centred at 2175 A (the UV bump, see Witt & 
Lillie |1973] and references therein) vanishes (almost) completely, 
while the slope of the extinction curve in the far-UV becomes 
steeper These differences depend on the mixture of dust grains 
and molecules present in a system. Tracing them through cosmic 
times improves our knowledge of the formation and evolution of 
dust. Furthermore, characterising the typical extinction curves 
of high-redshift galaxies is crucial for studies of the evolution of 
the SFR through cosmic times, in particular, if the SFR is deter- 
mined from the rest-frame UV stellar continuum. Therefore, the 
number of investigations on the properties of the dusty ISM at 
intermediate/high redshifts has been increasing in recent times 
(e.g., Malh otra.1997. Pitma n et al . |2000l MaioHno et al. 12001] 
Verne t et al. [ 20011 Vijh et al.l2003l Hopkins et al. 2004, Savaglio 
& Fall|2004l Wang et al.'2004', Wild & Hewett l2005.,NP05 l York 
et al. 2006). In particular, NP05 found that UV-luminous galax- 
ies at 2 < z < 2.5 can host dust-producing extinction curves with 
a significant UV bump. They made use of radiative transfer mod- 
els with the same dust/stars configuration that helps reproducing 
properties of dust attenuation in starburst galaxies in the local 
Universe (Gordon et al. 1997) as well as in the high-redshift one 
(Vijh et al. |2003l l. Under this assumption, NP05 concluded that 
the objects in their sample appear to host a mixture of dust grains 
that produce extinction curves analogous to those of the SMC 
and LMC and in between them. Interestingly, Wild & Hewett 
(I2005I I found that damped Lya absorbers with 0.84 < Zabs < 1 -3 
(identified through Ca 11 (H & K) absorption) exhibit extinction 
curves that are compatible with both the LMC and SMC ones. 



A possible connection between these two results is offered 
by our study. In fact, in the previous sections we have provided 
a robust evidence for the diversity of the extinction curves char- 
acterising 108 massive, UV-luminous galaxies at 1 < z < 2.5. 
Hence, the main conclusion of NP05 is extended down to z ~ 1 . 
This evidence stems from the distribution of the data points in 
the diagnostic diagram portraying the proxy for the strength of 
the 2175 A feature 734 versus the proxy for reddening in the rest- 
frame UV ySb (Fig. Is). In addition to the radiative transfer mod- 
els used by NP05 , we make use of additional ones that describe 
dust attenuation for a disc geometry and include the effect of an 
age-dependent extinction (Silva et al. 11998 ; Pierini et al. |2004| l. 
This analysis enforces the presence of a considerable fraction of 
dust out of the disc mid-plane of the sample galaxies, confirm- 
ing that the dust/stars configuration assumed by INP05I was not 
too wrong, though simple. The physical foundation of this com- 
plex dust/stars configuration is the existence of galactic winds 
(see, e.g., Murray et al. 2005). More noticeably, the evidence of 
dust with a UV bump also stems from the direct inspection of 
the composite spectra obtained for two subsamples of galaxies 
selected as a function of 734 either at 1 < z < 1.5 (Fig.|6]l or at 
2<z<2.5 (fig. 11 iniiPOS;. 



Our analysis clearly benefits from high S/N spec- 
troscopy, whereas others are limited by broad-band photome- 
try. Furthermore, unlike other samples, ours does not suffer from 
poor statistics and/or the removal of the dominant contribution 
to the emission by an AGN and/or the potential modification of 
the dusty ISM by the hard radiation field produced by AGN ac- 
tivity. Finally, it is made of galaxies detected in emission and not 
in absorption, which can thus be morphologically identified and 
studied at several wavelengths. These UV-luminous galaxies at 
1 < z < 2.5 are characterised by large total stellar masses and 



high SFRs as a consequence of the selection criteria (Sect. |2l0. 
Their stellar masses and SFRs tend to be lower the lower the 
redshift, consistently with the so-called "downsizing scenario" 
(Cowie et al. 1996, Gavazzi et al. 119961 Juneau et al. 2005ll. 
In fact, this scenario predicts a decrease with cosmic age of 
the transition mass separating actively from passively-evolving 
galaxies (e.g., Pannella et al. i2006i l. 

In our sample, galaxies that are heavily reddened in the rest- 
frame UV tend to be more massive and have higher SFRs, what- 
ever the redshift. Furthermore, they preferentially exhibit an ex- 
tinction curve with a significant UV bump. Nevertheless, we do 
not find any evidence that the observed strength of the UV bump 
is directly related to either stellar mass or SFR. Clearly these 
apparently conflicting results call for an interpretation that links 
properties of the stellar populations to those of the dusty ISM. 

For objects at 1 < z < 1.5 it is possible to measure the 
4000 A break index D4000 and/or the Balmer-break proxy C(36- 
39). These indices suggest that dust with a significant UV bump 
is present in galaxies where intermediate-age (i.e., from 0.2 to 
1 - 2Gyr-old) stellar populations exist and/or that are about 
three times older than those without a UV bump detection (see 
Sect. 14.4. 11 1. The existence of a link between extinction curve 
and the characteristic age of the stellar populations of a galaxy 
is nicely supported by the analysis of the absorption-line ratios 
that are available for the same objects. In fact, we find that the 
ratio of the absorption strengths in the Mgll doublet at 2800 A 
and in the Fe II doublet at 2600 A is larger in more UV-reddened 
galaxies with evidence of a UV bump. This holds whether the 
absorption arises in the ISM or in stellar atmospheres. Finally, 
galaxies with extinction curves exhibiting a significant UV bump 
at 1 < z < 1.5 predominantly appear as ordinary disc sys- 
tems, at variance with galaxies without a UV-bump detection at 
similar redshifts, which appear as irregular systems. Smoother 
morphologies are also consistent with a relatively longer star- 
formation history. 

In general, models with a mixed dust/stars configuration 
where MW-type dust is distributed in a geometrically thin disc 
(e.g., Silva et al. 119981 Pierini et al. 120041 ) can account for values 
of -2.5 < 734 < 1 in our sample of 1 < z < 1.5 galaxies (see 
Fig. |9l). However, they predict a rapid decrease of the observed 
strength of the UV bump as soon as the attenuation at UV wave- 
lengths becomes large along the line of sight. In order to describe 
the increase of the observed strength of the UV bump with in- 
creasing reddening and the existence of galaxies with 734 < -3, 
a screen-like geometry needs to be invoked (e.g., Witt & Gordon 
2000). The presence of large amounts of dust above the disc 
mid-plane can be obtained in principle via galactic winds, star- 
burst superwinds, and even AGN activity. It is favoured by the 
presence of magnetic fields (e.g., Greenberg et al. 119871 ). On the 
other hand, the ISM of objects with 734 ~ 1 does not seem to 
contain a significant fraction of carriers of the UV bump. The 
dust attenuation properties of these galaxies are similar to those 
of nearby starbursts (e.g., Calzetti et al. 119941 ). Consistently, our 
1 < z < 1.5 galaxies without a detection of the UV bump tend to 
be dominated by young (i.e., $ lOOMyr-old) stellar populations 
and to have irregular morphologies. 

There is no spectral indication for the characteristic ages of 
the stellar populations in the UV-bright galaxies at 2 < z < 2.5. 
However, the mean J -K colour of the subsample of FDF galax- 
ies at these redshifts ({J - K) - 1.7) suggests that the bulk of 



For instance, the need for spectra with a high S/N at rest-frame 
UV wavelengths implies a large fraction of young (i.e., < lOOMyr-old) 
stellar populations. 
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their stellar populations is most probably not much older than 
IGyr (Franx et al. "SOUS', van Dokkum et al. |2004l Pierini et al. 
I2005I I. Furthermore, their metallicities appear to be about solar, 
as well as for the galaxies at 1 < z < 1.5 (see Sect. |4.4.TT l. 

On the other hand, the subsample of UV-bright galaxies at 
2 < z < 2.5 offers the possibility to probe the ISM topology from 
the determination of the EWs of strong interstellar absorption 
lines (cf. Shapley et al. 2003). For these objects. .NP05 concluded 
that the presence of the UV bump does not seem to depend on 
the total metallicity, as given by the EW of the C IV doublet 
at 1550 A. Conversely, it seems to be associated with a large 
average EW of the six most prominent interstellar low-ionisation 
absorption lines in the far-UV. NP05 interpreted this result as 
indication of a larger covering fraction of young massive stars 
by neutral gas clouds (containing dust grains). 

At this point, it is possible to piece together the individual 
clues on the origin of the different extinction curves of UV- 
luminous galaxies at 1 < z < 2.5 and formulate a self-consistent 
scenario. In particular, we can address physical reasons for the 
presence or absence of the UV bump. First of all, we recall that 
the nature of the carriers of the 2175 A feature is probably mul- 
tiple, i.e., organic carbon (as in clusters of polycyclic aromatic 
hydrocarbons, PAHs) and amorphous silicates (Bradley et al. 
12005, and references therein). As supernovae (SNe) manufac- 
ture many of the condensible elements, clearly they are candi- 
date sources of dust. In particular, the expanding envelope of 
a type-II SN may represent an 0-rich environment where sili- 
cate dust can form. Silicon carbide and amorphous carbon are 
the main products of C-rich stars of intermediate mass (Whittet 
I2003I I. Conversely, intermediate-mass, O-rich stars produce sil- 
icate dust. These stars lose mass copiously during the red giant 
and asymptotic giant branch phases of their evolution. As seen 
before, AGB stars (with ages from 0.2 to 1 - 2 Gyr) are most 
probably present in large numbers in galaxies which exhibit the 
2175 A absorption feature in their spectra. These stars are con- 
sidered to be the main source of the probable carriers of the UV 
bump also in the local Universe (e.g., Galliano 2006jl. 

Once the basic units are produced and injected into the gen- 
eral ISM, dust reprocessing takes place there. Our study seems 
to confirm that dust self-shielding is important for the survival of 
the carriers of the UV bump in the general ISM (cf. Gordon et al. 
I2003D . despite other processes may modify or destroy them (see 
Whittet et al. l2004l l. This is consistent with recent studies on the 
formation and survival of PAH clusters in photodissociation re- 
gions (Rapacioli et al. 2006 and references therein). Milder UV 
radiation fields, as found in galaxies at 1 < z < 1.5, ease the 
need for self-shielding of course. An extreme example of this is 
the Milky Way, which has a relatively quiescent star-formation 
activity, a metal-enriched and aged ISM, and an average extinc- 
tion curve with a very pronounced UV bump (e.g., Cardelli et al. 
119891 see also Fitzpatrick & Massa l2007l for a recent analysis of 
the shapes of MW interstellar extinction curves). We recall that 
galaxies hosting dust with a significant UV bump at 1 < z < 1.5 
exhibit the heaviest UV reddening. In our preferred interpreta- 
tion, this is at least partly due to the screen-like effect of dust 
lifted from the disc mid-plane by galactic winds. 

Establishing the presence of different mixtures of dust grains 
in the dusty interstellar media of different high-redshift galaxies 
has important consequences on the perception of the cosmologi- 
cal evolution of the SFR. As shown in Fig. [15] attenuation mod- 
els where all stars are enshrouded by dust predict corrections of 
the luminosity at 1500 A (rest-frame) that differ by up to ~ 40% 
for the LMC and SMC extinction curves at fixed attenuation in 
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Fig. 15. Attenuation at 1500 A versus attenuation in the V band 
for three radiative transfer models of Witt & Gordon (2000) with 
clumpy SHELL dust geometry and the empirical Calzetti et al. 
(2000) law (see legend and Fig.[3]l. In the latter case filled circles 
are plotted in intervals of hE{B - V) = 0.1. 

the V-band (rest-frame) Ay < 2. Furthermore, a Calzetti law sys- 
tematically overpredicts the correction of L1500 and, thus, of the 
SFR with respect to an attenuation model with LMC-type dust. 
This is also true with respect to an attenuation model with SMC- 
type dust, the overestimate being much larger than 40% when 
Ay > 1 . In conclusion, not understanding amount, composition, 
topology, and three-dimensional distribution of dust in a galaxy 
can seriously undermine our knowledge of galaxy evolution. 

6. Conclusions 

High signal-to-noise optical spectra available for 108 massive, 
UV-luminous galaxies at 1 < z < 2.5 allow to constrain prop- 
erties of the extinction curve at rest-frame UV wavelengths as 
discussed by Noll & Pierini (2005. lNP05l l. As a main result, we 
find that these galaxies host different mixtures of dust grains and 
molecules, which produce different extinction curves. In partic- 
ular, the majority of the rest-frame UV spectral energy distri- 
butions (SEDs) suffering from heavy reddening exhibit a broad 
absorption excess centred at 2175 A (the so-called "UV bump"). 
Hence the associated galaxies host dust that can produce an ex- 
tinction curve similar to the average one determined for the in- 
terstellar medium (ISM) of the Large Magellanic Cloud or even 
of the Milky Way. Conversely, the extinction curve of galaxies 
with the least reddened UV SEDs is consistent with the average, 
featureless, and steep extinction curve of the Small Magellanic 
Cloud, if the strength of the apparent UV bump is not consider- 
ably reduced by dust-stars mixing and age-dependent extinction 
effects as suggested for local disc galaxies. 

For objects at 1 < z < 1.5, the comparison of measured 
and synthetic stellar spectral indices shows that dust with a sig- 
nificant UV bump is present in galaxies that host a rather large 
fraction of intermediate-age stars (i.e., from 0.2 to 1 - 2 Gyr old) 
at these redshifts. The older ages and the higher reddening at 
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UV wavelengths found for these galaxies are consistent with the 
larger Mg II-to-Fe II ratio that we measure. At the same time, 
older ages are consistent with a larger fraction of smooth, radial 
surface brightness profiles in the rest-frame optical. Therefore, 
at 1 < z < 1.5 galaxies with evidence of a UV bump appear 
as disc systems which are more evolved than the more compact, 
clumpy objects whose dust properties resemble those of nearby 
starburst galaxies. Nevertheless, part of the dusty ISM seems to 
be lifted above the disc mid-plane owing to the action of galactic 
winds. 

In a complementary way, at least for objects at 2 < z < 2.5, 
the equivalent width of the most prominent interstellar absorp- 
tion lines in the far-UV carries information on the topology of 
the ISM. As INP051 showed, the covering fraction of young mas- 
sive stars by neutral gas clouds (containing dust) is larger in 
galaxies with evidence of a UV bump than in those with an un- 
detected 2175 A feature. 

The previous results can be interpreted in a self-consistent 
way by a scenario where the basic units of the (probably multi- 
ple) carriers of the UV bump mostly originate from intermediate- 
mass (2-5 Mq) stars in late evolution phases. A high production 
rate of the carriers is, therefore, connected to high amounts of the 
required stars, which needs suitable star-formation histories. The 
carriers of the 2175 A feature are most probably made of organic 
carbon and amorphous silicates. They can survive in the harsh 
environments of UV-luminous high-redshift galaxies owing to 
dust self-shielding, which is most efficient in massive, gas-rich, 
high-metallicity galaxies. As indicated by our analysis, good 
candidates for galaxies with a pronounced UV bump in their ex- 
tinction curves (obviously fulfilling the requirements above) are 
highly dust-reddened, UV-bright galaxies at 2 < z < 2.5 and 
near-IR bright, intensely star-forming galaxies at 1 < z < 1.5. 

The existence of different extinction curves implies that dif- 
ferent patterns of evolution and reprocessing of dust exist among 
massive, UV-luminous galaxies at high redshift. Ignoring this 
may produce a non-negligible uncertainty on the SFR estimate 
of a galaxy based on the rest-frame UV. 
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Table 1. Basic parameters, i.e., ID, world coordinates (J2000), R and A's magnitudes, and redshift, of the 108 sample galaxies from 
the FDF, K20, and GDDS. For the origin of the data see Sect. |2] Furthermore, the main parameters of this study, j3b and 734 (see 
Sect. 13. Il l, are listed. Due to non-optimal wavelength ranges, fii, could only be measured for 88 objects. 
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Table 2. Average properties of the six subsamples defined in Sect. 14.41 The numbers of objects used for the derivation of the 
morphological parameters (in parentheses) are smaller than for the complete subsamples, since Hubble ACS images are not available 
for the K20 quasar field (Sect. l2!2] i and some outer parts of the FDF (Sect. lZTT i. In the following we list the magnitudes/? and K^, the 
UV-continuum slope parameter y6b, the UV-bump tracer 734, the luminosity at 1500 A Ligoo, the dust reddening parameter E(B - V), 
the bolometric luminosity Lboi and the SFR from the model fits, the total stellar mass Mstars, the specific SFR <p, the effective radius 
Rs, the Sersic index «ser, the ellipticity e, the CAS parameters concentration C and asymmetry A, and the "dumpiness" R-y + Ra- 



Parameter 


1 < z < 1.5 


1 < z < 1.5 


1 <z 


< 1.5 


2 < z < 2.5 


2 < z < 2.5 


2 < z < 2.5 




Pb < -1-5 


ySb > -1-5 


j8b > 


-1.5 


j8 < -0.4 


j8 > -0.4 


P > -0.4 






734 > -2 


734 < -2 




734 > -2 


734 < -2 




13(12) 


14(12) 


13(9) 


16 


9 


9 


R Tniasl 


23.44 ±0.11 


23.17 ± 0.09 


23.21 


±0.17 


23.89 ± 0.15 


23.77 ± 0.13 


23.69 ± 0.13 


rmasl 


20.94 ± 0.24 


20.24 ±0.11 


19.73 


±0.09 


21.22 ±0.18 


20.66 ± 0.22 


20.52 ± 0.20 




-1.70 ±0.06 


-0.95 ±0.11 


-0.72 


±0.10 


-1.96 ±0.09 


-1.61 ±0.16 


-1.37 ±0.15 


734 


-0.44 ±0.18 


-0.63 ± 0.29 


-3.70 


±0.31 


+0.36 ±0.36 


-0.44 ± 0.35 


-3.50 ±0.30 


logLisoo [W/A] 


33.69 ± 0.07 


33.67 ± 0.06 


33.48 


±0.07 


34.10 ±0.05 


34.06 ± 0.06 


34.17 ±0.05 


Eb-w 


0.26 ±0.01 


0.38 ± 0.02 


0.42: 


t0.02 


0.27 ± 0.01 


0.37 ± 0.02 


0.39 ± 0.02 


logLbo, [W/A] 


11.55 ±0.09 


12.02 ±0.06 


11.99 


±0.07 


11.97 ±0.08 


12.39 ±0.08 


12.56 ±0.08 


logSFR[Mo/yr] 


1.83 ±0.08 


2.30 ± 0.08 


2.23 : 


t0.09 


2.24 ± 0.09 


2.61 ±0.08 


2.82 ± 0.08 


logM,,,,, [Mq] 


9.91 ±0.21 


10.18 ±0.09 


10.40 


±0.06 


10.14±0.11 


10.67 ±0.11 


10.76 ±0.13 


log0[Gyr-'] 


0.89 ±0.17 


1.11 ±0.14 


0.79: 


t0.16 


1.06 ±0.13 


0.93 ±0.11 


1.06 ±0.14 


i?c [kpc] (from log i?e) 


2.85 ± 0.43 


3.59 ±0.33 


3.99: 


tO.52 


2.79 ±0.52 


3.73 ± 1.36 


3.83 ±0.82 




1.32 ±0.36 


0.98 ± 0.23 


0.54: 


tO.25 


2.07 ± 0.40 


2.48 ± 0.47 


1.61 ±0.49 


e 


0.49 ± 0.07 


0.42 ± 0.06 


0.45 : 


t0.07 


0.63 ± 0.04 


0.49 ± 0.07 


0.62 ± 0.07 


CcAS 


2.43 ± 0.09 


2.16±0.11 


2.07 : 


tO.U 


2.52 ±0.11 


2.46 ±0.12 


2.19 ±0.17 


^CAS 


0.20 ± 0.06 


0.19 ±0.04 


0.28: 


t0.08 


0.18 ±0.02 


0.18 ±0.05 


0.22 ± 0.05 


Rj + Ra 


0.24 ± 0.02 


0.24 ± 0.03 


0.26: 


t0.04 


0.25 ± 0.05 


0.40 ± 0.09 


0.32 ± 0.05 



